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Abstract —The Delay Constrained Relay Node Placement 
(DCRNP) problem in Wireless Sensor Networks (WSNs) aims 
to deploy minimum relay nodes such that for each sensor node 
there is a path connecting this sensor node to the sink without 
violating delay constraint. As WSNs are gradually employed 
in time-critical applications, the Importance of the DCRNP 
problem becomes noticeable. For the NP-hard nature of DCRNP 
problem, an approximation algorithm-Set-Covering-based Relay 
Node Placement (SCA) is proposed to solve the DCRNP problem 
in this paper. The proposed SCA algorithm deploys relay nodes 
iteratively from sink to the given sensor nodes in hops, i.e., in the 
fcth iteration SCA deploys relay nodes at the locations that are 
k hops apart from the sink. Specifically, in each iteration, SCA 
first finds the candidate deployment locations located within 1 
hop to the relay nodes and sensor nodes, which have already 
been connected to the sink. Then, a subset of these candidate 
deployment locations, which can guarantee the existence of paths 
connecting unconnected sensor nodes to the sink within delay 
constraint, is selected to deploy relay nodes based on the set 
covering method. As the iteration of SCA algorithm, the sensor 
nodes are gradually connected to the sink with satisfying delay 
constraint. 

The elaborated analysis of the approximation ratio of SCA 
algorithm is given out, and we also prove that the SCA is a poly¬ 
nomial time algorithm through rigorous time complexity analysis. 
To evaluate the performance of the proposed SCA algorithm, 
extensive simulations are implemented, and the simulation results 
show that the SCA algorithm can significantly save the deployed 
relay nodes comparing to the existing algorithms, i.e., at most 
31.48% deployed relay nodes can be saved due to SCA algorithm. 

Index Terms —Wireless sensor networks, set covering, delay 
constrained relay node placement, approximation algorithm. 

I. Introduction 

IRELESS Sensor Networks (WSNs) attract consid¬ 
erable attention in recent years for their immense 
potential in the applications of environment monitoring, battle- 
held reconnaissance and industrial automation ID-El. WSNs 
comprise spatially distributed Sensor Nodes (SNs), which 
sense certain circumstance information, and one or several 
sinks to collect the information sensed by SNs. Typically, SNs 
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are cheap and powered by batteries. This results in that the 
power and communication radii of a SN is limited. Thus, to 
prolong the network lifetime and improve the scalability, Relay 
Nodes (RN) are introduced into WSNs. RNs are equipped 
with plenty power and large communication radii, however, the 
RNs are also costly Ei-il. Hence, the Relay Node Placement 
(RNP) problem try to build connected WSNs by deploying a 
minimum number of RNs subject to various constraints, such 
as network lifetime, throughput and delay. 

Recently, the importance of Delay Constrained RNP 
(DCRNP) problem is highlighted by the employment of WSNs 
in the time-critical applications, e.g., the factory automation 
and patent monitoring S-El. In the factory environment, the 
data sensed by SNs is typically time-sensitive, such as alarm 
notification and information for feedback control, and thus the 
importance of receiving the data at the sink in a timely manner 
is noticeable Q-llll. It has been reported in ifT^ that DCRNP 
problem is NP-hard, which implies that no polynomial time 
algorithm exists for this problem as long as P ^ NP. 

Although RNP problem has been extensively studied, the 
literature about DCRNP problem is rare. Bhattacharya and 
Kumar lIT^ - ll^ first proved that DCRNP problem is NP- 
hard, and then proposed a polynomial time approximation 
algorithm to solve this problem. Nigam and Agarwal ll23ll 
designed an algorithm to optimally solve DCRNP problem 
without guaranteeing a polynomial time complexity. As in 
ca-Ei, in this paper, RNs can only be placed at a subset of 
the predetermined locations, which are called the Candidate 
Deployment Locations (CDLs). 

In order to devise a polynomial time algorithm with better 
performance, a Set-Covering-based Approximation (SCA) al¬ 
gorithm is proposed in this paper to yield a desirable solution 
for DCRNP problem. The proposed SCA is an iterative 
algorithm and deploys RNs in hops from the sink to SNs. 
In order to maintain network connectivity, in each iteration, 
the RNs that are deployed based on set covering method 
have 1-hop neighbor SNs or RNs, which have been connected 
to the sink within delay constraint. As SCA algorithm is 
iteratively executed, SNs are gradually connected to the sink 
with fulfilling delay constraint. To evaluate the proposed SCA 
algorithm, extensive simulations are implemented, and the 
simulation results show that the SCA algorithm outperforms 
the existing algorithms outstandingly. 

The major contributions of this paper are listed as follows: 

• An approach is designed to represent each CDL or SN by 
the unconnected SNs, which can be connected to the sink 
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by the paths that fulfil delay constraint and pass through 
this CDL or SN. 

• An approximation algorithm-SCA is proposed for 
DCRNP problem. In each iteration, SCA selects a subset 
of the CDLs, which are located within 1-hop from the 
RNs and SNs that have been connected to the sink, to 
deploy RNs by solving a set covering problem that are 
modeled by the approach mentioned above. Due to this 
manner, SCA can avoid the limitation suffering by the 
algorithm presented in ifTOl - OTll . 

• The approximation ratio of SCA algorithm is explicitly 
analyzed, and the time complexity of SCA algorithm is 
also detailed in this paper. 

This paper is organized as follows. Section II reviews 
the related works. The problem formulation is given out in 
Section III. The preliminaries and heuristic are given out in 
Section IV. Section V describes our proposed algorithm for 
the DCRNP problem. Section VI shows an explicit analysis 
of approximation ratio and time complexity of the proposed 
heuristic algorithm. The efficiency of the proposed algorithm is 
validated through extensive simulations in Section VII. Finally, 
Section VIII concludes the whole paper. 

II. Related Work 

The related literatures are divided into two categories, 
i.e., the literatures about the DCRNP problem and the RNP 
problem without delay constraint. 

A. RNP Problem without Delay Constraint 

Lin and Xue ID studied the RNP problem and formulated 
it as the as the Steiner minimum tree with minimum number 
of Steiner points and bounded edge length problem (SMT- 
MSP). Then, they proved the SMT-MSP to be NP-Complete 
and proposed a 5-approximation algorithm. Chen et al. Col 
demonstrated that the algorithm proposed in 13 is actual a 
4-approximation algorithm, and proposed a 3-approximation 
algorithm for the RNP problem. Cheng et al. im presented 
a 3-approximation algorithm and a 2.5-approximation algo¬ 
rithm based on so-called 3-star structure. Lloyd and Xue 
II3 presented a 7-approximation algorithm for the single- 
tired network RNP problem and a (5 -I- e)-approximation 
algorithm for the two-tired network RNP problem. Srinivas et 
al. m studied the problem of constructing and maintaining 
the wireless backbone network for the WSNs, and they also 
considered mobile wireless backbone network. 

Misra et al. IfT^ - lfTTl studied the constrained RNP problem 
in the single-tiered network and proposed a polynomial time 
0(l)-approximation algorithm. Yang et al. d first stud¬ 
ied the constrained RNP problem in the two-tiered network 
and proposed an algorithm with 0(l)-approximation ratio 
for 1-connected single cover problem. Then they proposed 
two algorithms with 0(l)-approximation ratio and O(lnn)- 
approximation ratio under different settings for 2-connected 
double cover problem. 


B. DCRNP Problem 

The DCRNP problem has been studied in lfT9l - ll2^ . Bhat- 
tacharya and Kumar lfT9l - ll2Tll first proved that the DCRNP 
problem is NP-hard, and a shorted path tree based algorithm 
is proposed. This algorithm preliminarily forms a shortest 
path tree to connect each SN to the sink, then it saves 
the deployed RNs by gradually removing the RNs on the 
shortest path tree. This brings in a limitation that the deployed 
RNs can only be those contained in the originally yielded 
shortest path tree, and the worst case of this algorithm happens 
when all the RNs of the optimal solution have been missed 
by the shortest path tree. Sitanayah et al. i23i studied the 
fault-tolerant RNP problem under the length constraint, and 
proposed a heuristic to solve this problem. However, no time 
complexity analysis and performance guarantee were given. 
Nigam and Agarwal 12^ formulated the DCRNP problem as 
a linear programming problem, and proposed a branch-and-cut 
algorithm to optimally solve the DCRNP problem. However, 
the proposed algorithm can only solve a special case of the 
DCRNP problem (each of the source node cannot have a 
singleton node cut), and the time complexity of this algorithm 
grows exponentially, which indicates it cannot be applied to 
the large scale problem. 

Some researches on the related scopes are also listed in 
follows. Voss II 24 I used the hop count to represent the network 
delay, and studied the hop constrained Steiner tree problem. 
A local search algorithm based on tabu search was proposed, 
but this paper did not analyze the approximation ratio and 
the time complexity. Costa et al. 1251 studied the Steiner tree 
problem with revenue, budget and hop constraints. A tabu 
search based heuristic was proposed to solve this problem, 
but the proposed heuristic is not a polynomial time algorithm. 
Gouveia et al. 12^ studied the distance constrained minimum 
spanning tree problem, and three approaches were proposed 
based on the column generation scheme and the Lagrangian 
relaxation. However, no time complexity analysis is provided 
in the paper. 

III. Problem Formulation 

The network delay is typically composed of the processing 
delay, the queuing delay, the transmission delay and the 
propagation delay. In the WSNs, the distance between two 
different SNs is normally a few tens or hundreds of meters. 
Therefore, normally compared to the other three kinds of 
delay, the transmission delay can be ignored. Since the rest 
of delays are at proportional to the hop count, in this paper, 
the network delay is represented by the hop count. Therefore, 
DCRNP problem is reduced to a hop constrained RNP problem 
in the rest of this paper. 

Given a set of SNs S = {si, S 2 ) s„}, a set of Candidate 

Deployment Locations (CDLs) C = {ci, C 2 ,..., Cm} for de¬ 
ploying RNs and a sink K, we can build an undirected graph 
G = {V,E}, where V = -Sy CyiiT} is the node set, and 
E is the edge set. Vu, v € V (u ^ v), if u and v are the 
two ends of an edge in E, u and v should fulhl the following 
conditions: 
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• If u S S' or z; G S', then u and v should meet that ||u — 
z;|| < r; 

• if u ^ S and v ^ S, then u and v should meet that 

||u — ull < R, 

where ||m — u|| denotes the Euclidean distance between u and 
V. r and R (r < R) are the communication radii of the SN 
and RN, respectively. Besides, in this paper a path between u 
and V is denoted by p{u, v). 

Definition 1 (DCRNP Problem): For an undirected graph 
G = {V,E}, the DCRNP problem searches for an induced 
subgraph G' = {V',E'} of G, where V = SIJG'IJI^}, 
and G' is a subset of C with the minimum cardinality such 
that the following condition is satisfied: there exists at least 
one path, which complies with the delay constraint, between 
each SN and the sink. 

Let Pi = {pi{si,K),p 2 {si,K), ...,pufisi,K)] (1 < i < n) 

be the ki paths connecting the SN Si and the sink K, where 
Pj{si,K) (1 < j < ki) denotes the jth path between Si and 

n 

K. Let P = yj Pi denote the set of all the paths connecting 

2—1 

the given SNs and the sink K. 

CDLs on the induced subgraph G' are selected to deploy 
RNs, thus, the terms RN and CDL are interchangeable in this 
paper. Consequently, a node on G' can be a SN, RN or the 
sink. Let N{p) denote all the nodes on path p. Then, let C(p) 
denote RNs on path p, and the notation C{p) represent the hop 
count of path p. Obviously, we have that C{p) = |A/"(p)| — 1. 

The DCRNP problem can be formulated as an optimization 
problem: 

Minimize \G'\ (la) 

s.t. VP* G P (1 < i < n), 

3pj{si,K) G Pi (1 < j < \Pi\),C{Pj{si,K)) < Ai 

where Ai is the delay constraint for the *th SN. Without loss 
of generality, in this paper we assume that the given SNs have 
the same delay constraint, i.e., Ai = A, z G {1, 2,..., n}. 

The DCRNP problem has been proved NP-hard IT^ , there¬ 
fore, this paper proposes a heuristic algorithm to approxi¬ 
mately solve the DCRNP problem. The path connecting the 
zth SN Si and the sink and satisfying the delay constraint A 
is called a feasible path of Si, and the subgraph satisfying 
the formulation (fTbl i is called a feasible graph of the DCRNP 
problem. 

IV. Preliminaries and Heuristic 

Feasible graphs for a DCRNP problem should be connected. 
Thus, we can find an induced subgraph, which is a tree rooted 
at the sink and connecting all the SNs, from each feasible 
graph. Therefore, in this paper, the proposed algorithm yields 
a feasible graph that is a tree for each given DCRNP problem. 

Given a tree T and two different nodes of T, u and v, 
we denote pt{u,v) as a path from m to u of T. Vs G S, 
if C{pt{K, s)) < A, we call T a feasible tree for the given 
DCRNP problem. 

Next, the conception of the level of a tree will be introduced 
to facilitate the explanation of the proposed SCA algorithm. 
Given a feasible tree T, the level of a node g in T is the 



Fig. 1. An illustration of the conception of the level. 


hop count of the path between the sink K and node q, 
i.e., C{pT{K,q)). The set of nodes on the fcth level of T 
is denoted by L^. As illustrated in Fig. [T] = {cijCs}, 

Lt = {c2,C3,ce}, = {si,S 4 ,C 4 ,C 7 } and = { 52 , 53 }. 

Observing that the nodes on each level should have 1-hop 
neighbors on their adjacent levels to maintain the network 
connectivity, we can deploy RNs from sink to SNs level by 
level such that the sink is gradually connected to SNs by the 
deployed RNs. Suppose that we have deployed RNs from the 
1st level to the (fc — l)th level of T, and now attempt to deploy 
RNs for the fcth level. Let denote the set of 1-hop 

neighbors of the nodes in Obviously, the nodes on fcth 

level should be 1 -hop neighbors of nodes on (fc — l)th level, 
i.e., C 

Therefore, we move on to search the 1-hop neighbors of the 
nodes on (fc — l)th level. As the search is accomplished, i.e., 
is known, we try to select the nodes on fcth level 
from V(Ly~^). The nodes that are selected to be located on 
the fcth level should guarantee the existence of feasible paths 
for the unconnected SNs. For this purpose, an approach, which 
represents each node q in by the unconnected SNs 

that can be connected to the sink by the feasible paths passing 
through q, is designed and shown in the following path of this 
section. 

In this paper, the shortest path between u and v means a 
path connecting u and v with least hop count, and the notation 
'H{u,v) stands for a shortest path between u and v. Let Pi 
(1 < i < n) be the set of all feasible paths between a 
SN Si and the K, i.e., \IC{p{si,K)) < A, p{si,K) G Pi. 
Consequently, the set of nodes lying on the feasible paths of 
Si is y N [x). Let q be an arbitrary CDL or SN. 

X^Pi 

Theorem 1: The sufficient and necessary condition for q G 
y Af (x) is that 

xePi 

C{-H{q,Si))pC['H{q,K))<A. (2) 

Proof: Sufficiency: From g G y Af (a;) we knows that 

X^Pi 

there is at least one path, p (p € Pfi, passing through q, 
therefore, we have C{p) < A. 

Let p (p) be a segment of p between q and Si (K). Then, 
we have 


C{p) > C{'H{q,Si)) 


(3) 
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and 

C{p)>C{H{q,K)). (4) 

Since p consists of p and p, p is a feasible path, which 
implies 

C(p) = C(p)+C(p) < A. (5) 

Plunging inequalities Q and (IDl into (|5]l, we can achieve 

C{'H{q,Si))+C{n{,q,K))<l^. ( 6 ) 

So far, we have completed the first part of the proof. 

Necessity; We build two shortest paths which connect q 
and Si and connect q and K, respectively, i.e., H(g, Si) and 
'H{q,K). Then, we can form a path p between Si and K by 
combining (g, Si) and 7 ^( 9 , K). According to the inequality 
0 , the hop count of p is given by 

C{p) = C{n{q,s{))+C{n{q,K)) 

< A, ^ 

which confirms that p is a feasible path between Si and K. 
Therefore, we can conclude that 

g e IJ A/'(x). (8) 

X^Pi 

This completes the proof of Theorem [T] ■ 

Theorem [T] implies that if a RN or SN q lies on a feasible 
path of a SN, then q should satisfy inequality 0. Let Sk be 
the set of unconnected SNs, which are located on the levels 
higher than k. As mentioned above, we expect to select the 
nodes on the fcth level of T from such that feasible 

paths exist for the unconnected SNs. To this end, for each node 
q in V(Ly“^), we try to find the unconnected SNs, which can 
be connected to sink by the feasible paths passing through q, 
and the set of these SNs is denoted by Q{q). According to 
Theorem [T] the SNs in Q{q) should fulfil that 

C(7f(s,g))+fc< A, (9) 

where s G Q(q), and we say that s is connected to the sink 
by q. 

When q (q G is selected as a node on the fcth 

level, according to Theorem [T] we can promise the existence 
of feasible paths, which pass through q and connect the SNs 
in Q{q) to the sink within delay constraint. 

However, the inequality 0 cannot ensure that the RNs 
deployed on current level are closer to the unconnected SNs 
than the nodes on previous levels, which may bring in a large 
amount of redundant deployed RNs. As shown in Fig.|2] delay 
constraint is set as A = 7. According to inequality 0, we can 
obtain that Q(ci) = {si, 32,33} and Q(c2) = {32,33}. Since 
V(L!}) = {ci, C2} and all the SNs can be connected to the sink 
by Cl, to select as fewer CDLs as possible on each level, only 
Cl is selected to be allocated on the 1st level. This ignores that 
the feasible path that passes through ci and connects S2 and 
33 to the sink is much longer than the feasible path passing 
through C 2 to S 2 and 33 , i.e., C(p(ci, 32 )) > C(p(c 2 , 32 )) and 
C(p(ci,S 3 )) > C(p(c 2 ,S 3 )). This ignorance leads to a local 
optimum, whose edges are denoted by the blue dashed line 
segments in Fig. |2] and misses the optimal solution, whose 



Fig. 2. An illustration of the redundant RNs brought by the utilization of 
inequality 


edges are represented by the black line segments. As a result, 
totally 3 redundant RNs are introduced. 

To this end, we introduce a rule to guarantee that the 
RNs deployed on each level gradually move closer to the 
unconnected SNs. Let g be a node on the (fc — l)th level, 
and g be a 1 -hop neighbor of g. 

Lemma 1: Q{q) C Q{q). 

Proof: For each SN in Q(g), the inequality 0 is met, 

i.e., 

V 3 G Q(g), C(7f(s, g)) + fc < A. (10) 

Since g is a 1-hop neighbor of g and g is a node on the (/c—l)th 
level, the following inequality is fulfilled 

V3 G Q(g), 3 p( 3 ,g), C(p(s, g))-b fc - 1 < A. (11) 

Due to the fact that C(p( 3 , g)) > C{'H{s, g)), we can conclude 
that 

VsGQiq), C{n{s,q)) + k-l<A. (12) 

which implies that V3 G Q(g), s G Q(g). Thus Q(g) is a 
subset of Q(g). This completes the proof. ■ 

According to Lemma [T] the SNs in Q(g) can also be 
connected to the sink by g. Therefore, to guarantee that the 
nodes on the kth level are closer to the unconnected SNs than 
their 1 -hop neighbors on the (fc — l)th level, we restrict the 
nodes on the fcth level to satisfy a rule 

C( 7 f(g, 3 )) <C( 7 f(g, 3 )), (13) 

where 3 G Q{q) and s G Q(g). Finally, the SNs in Q{q) 

should meet both the inequalities 0 and (foT l. 

As illustrated in Fig. |3l ci is the only node on the 1st 
level, and 3 i is connected to the sink by ci. Additionally, 
the delay constraint is set as 6 , and C('H(ci, 3 i)) = 4. The 
1 -hop neighbors of ci are represented by the green circles, 
i.e., C 2 , C 3 , C 4 , C 5 and cq. The shortest paths from these 1- 
hop neighbors to 3 i are represented by the green dashed line 
segments. It can be seen from Fig. |3that the inequality 0 is 
met by these 1 -hop neighbors except for ce, but only C 3 and 
C 4 have a hop count to si less than 4. Thus, in these 1-hop 
neighbors, only C 3 and C 4 fulfil both the formulations 0 and 
(HI, i.e., Si can be connected to the sink only by C 3 and C 4 . 

To connect all the unconnected SNs to the sink, we should 
select a subset U of such that each unconnected SN 

can be connected to the sink by these nodes, i.e., (J Q{u) = 

ueu 
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Fig. 3. An illustration of how to determine the SNs connected to the sink by 
a CDL. 


Sk- Furthermore, to save the deployed RNs we expect to min¬ 
imize the cardinality of subset U. Thus, we can formulate the 
problem, which requires to select the nodes on kth level from 
V(Ly~^), as the set covering problem described follows; given 
an unconnected SNs set Sk and a node set whose 

each element q is represented by the unconnected SNs that can 
be connected to the sink by q, then this problem searches for 
a minimum subset U of such that [J Q{u) = Sk- 

u&U 

As a result, we can employ the set covering algorithm to solve 
this problem. Due to the effective performance and the wide 
utilization, the greedy heuristic algorithm is employed to solve 
this set covering problem. 

V. Description of SCA 

The proposed SCA algorithm is composed of three steps and 
detailed in Algorithm [T| In the hrst step, the SCA algorithm 
preliminarily tests either a feasible solution exists or RNs are 
necessary for the DCRNP problem based on shortest path tree 
algorithm. 

The second step is the principal body of the SCA algorithm. 
The second step of SCA is iteratively executed so as to allocate 
the nodes on each level. To be specihc, in the fcth iteration, 
SSCA moves on to select nodes on the fcth level. Firstly, 
the second step of SCA searches for 1-hop neighbors of the 
nodes on the {k — l)th level, i.e., Next, for each 

node q in V(Ty~^), the second step of SCA searches for 
the unconnected SNs, which can be connected to the sink by 
q, i.e., Q{q). Then, each node q in is represented 

by Q{q), and as a result the current problem is transformed 
into the set covering problem, which attempts to select a 
minimum subset of to fully cover the unconnected 

SNs. Finally, the greedy heuristic algorithm is employed to 
solve this set covering problem. Consequently, the selected 
subset of is the set of nodes on fcth level. As the 

execution of the second step of SCA, a feasible tree T is 
gradually built, and clearly, due to the delay constraint, the 
number of iteration of the second step of SCA cannot be larger 
than A. 

The process of the second step of SCA is illustrated in Fig. 
m where the circles and triangles with the same color are the 1 - 
hop neighbors of the nodes on the same level, and the circles 


Algorithm 1 SCA. 

Require: 

A set of SNs S = {si, S 2 ? •••) a set of CDLs C = {ci, C 2 ,Cm}, 
a sink node K, the communication radii of the SN r and the RN R, the 
delay constraint A. 

Ensure: 

A feasible tree T. 

1: trap =a shortest path tree rooted at K and including all SNs in S with 
using CDLs in C; %the first step of SCA begins. 

2: if Vi G {1, 2, ptmp{si,K) < A then 

3: tmp =a shortest path tree rooted at K and including all SNs in S 

without any CDLs in C; 

4: if 3i G {1,2, pt 7 np(si,K) > A then 

5: input S, C and K into the second step of SCA; %the first step of 

SCA ends. 

6: C' =a subset of C returned by the second step of SCA; 

7: input S, C and K into the third step of SCA; 

8: T =a feasible tree returned by the third step of SCA; 

9: else 

10: declare that no RNs are necessary to build a feasible tree T and 

terminate SCA; 

11: end if 

12: else 

13: declare that no feasible solutions exist for this DCRNP problem and 

terminate SCA; 

14: end if 
15: return T; 



Fig. 4. An illustration of the procedure of the second step of SCA. 


and triangles with red edges are the SNs or RNs on each 
level. At the beginning, the second step of SCA first searches 
the 1 -hop neighbors of the sink K, and these neighbors are 
denoted by the black circles and triangles. Next, the second 
step of SCA represents these neighbors by the unconnected 
SNs, and the greedy heuristic is employed to hnd a set cover 
for the unconnected SNs. As shown in Fig.|4l c\ is selected in 
this iteration. In the second iteration, the second step of SCA 
hrst hnds the 1 -hop neighbors of ci, and these neighbors are 
denoted by the blue circles and triangles. Then, s, and C 3 
are selected in this iteration. As the iteration repeats, all the 
unconnected SNs are gradually connected to the sink. The 
second step of SCA is detailed in the Algorithm |2] 

Remark 1: As 1-hop neighbors have been represented by 
the SNs connected to the sink by them, the situation that 
several neighbors connect the same unconnected SNs to the 
sink may happen, i.e., Q{q) — Q{q), where q and q are all 
1-hop neighbors of the nodes on the last level. To deal with 
this situation, the weight of a neighbor, q, is introduced and 
defined as follows. Let Q{q) be the set of unconnected SNs, 
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which are connected to the sink by q. Then the weight, uj{q), 
of q is given by 


U}{q) 


U C('W(g,x)) 

x&Q{q) 


(14) 


If the situation mentioned above happens, we select the one 
with least weight. 


Algorithm 2 The Second step of SCA. 

Require: 

A set of SNs S = {si, S 2 i •••i Sn}, a set of CDLs C = {ci, C 2 , 
a sink node K, the communication radii of the SN r and the RN R, the 
delay constraint A. 

Ensure: 

A set, C', of CDLs. 

1: V(Ly) = the 1 -hop neighbors of JC in S and C; 

2: LO, = K- 
3: fc = 1; 

4: Sk = S; 

5: while (5 ^ 0) do 
6: Sfc = Sfe-V(LO); 

7: i = 0; 

8 : tmp =the 2 th element in V(-L^); 

9: while {tmp ^ 0) do 

10: con[i] = Q{tmp)', 

11 : weight[i] = <jj{tmp)', 

12 : 2 = 2 + 1 ; 

13: tmp =the 2 th element in V(-L^); 

14: end while 

15: sort V(L^) according to the cardinality of each cor 2 [ 2 ] and the 

weight[i]', 

16: = a subset of V(-L^) found by the greedy heuristic algorithm; 

17: Sk = Sk- 

18: C = C-L%- 

19: V(Ly) = the 1-hop neighbors, which are selected from Sfc and C, of 

the nodes in 

20: tmpRe = — S', 

21 : C = C \JtmpRe', 

22: fc = fc -f 1; 

23: end while 
24: return C"; 


The redundant RNs may be deployed by the second step 
SCA for the introduction of the rule (fOl) . which is illustrated 
in Fig. |5] The feasible tree returned by the second step 
SCA is denoted by the black line segments in Fig. |5] Since 
C{'H{K, si)) < C{'H{c 3 , si)), Si cannot be connected to the 
sink by C 3 , and the second step SCA select two CDLs to 
deploy RNs, i.e., ci and C3, on the first level. However, there 
is a feasible path passing through C 3 and connecting K to si. 
Thus, the redundant RNs Ci and C 2 are introduced. Actually, 
a solution with fewer RNs exists in Fig. and this solution 
is represented by the red dashed line segments. 

To save the redundant RNs introduced by the second step of 
SCA, the third step of the SCA algorithm is designed. Since 
we expect that the number of RNs deployed on each level is 
as fewer as possible, we try remove the RNs connecting fewer 
SNs to the sink on each level. Thus, each RN q returned by the 
second step of SCA is assigned a weight, which is defined as 
the number of SNs meeting the inequality (|2|l for q. Then, the 
returned RNs are sorted in ascending order according to the 
weight, and we try to gradually delete the redundant RNs from 
the one with least weight. As a RN is deleted, the shortest path 
tree is formed to check whether a feasible tree exists for the 



Fig. 5. An illustration of the redundant CDLs introduced by the second step 
SCA. 


DCRNP problem on the remaining RNs. If such feasible tree 
exists, the RN is deleted forever, otherwise the RN is retrieved 
and marked as tried. This procedure is repeated until all the 
RNs is deleted or marked as tried. The third step of the SCA 
algorithm is detailed in the Algorithm [3] 


Algorithm 3 The Third step of SCA. 

Require: 

The set S of SNs, a set C' of CDLs returned by the second step of SCA. 

Ensure: 

A feasible tree T. 

1: T = a shortest path tree rooted at the sink and including all the given 
SNs by using the CDLs in C'; 

2: for all tmp £ C do 
3: calculate the weight of imp; 

4: end for 

5: sort the CDLs in C in an ascending order according the weight; 

6 : tmp = the first element of C'', 

7: C' = C' - tmp; 

8 : while {tmp 5 ^ 0 ) do 

9: tmpT = a shortest path tree rooted at the sink and including all the 

given SNs by using the CDLs in C'; 

10: if {tmpT is a feasible tree) then 

11: T = tmpT; 

12: C = the CDLs on T; 

13: sort the CDLs in C' in an ascending order according the weight; 

14: tmp = the first element of C'; 

15: C = C — tmp; 

16: else 

17: mark tmp as tried; 

18: C = C'\Jtmp; 

19: tmp = the least weighted untried CDL in C'; 

20: end if 

21: end while 
22: return T; 


VI. Analysis of SCA 


A. Approximation Ratio 

Let T* be a optimal feasible tree for the given DCRNP 
problem, and the set of the RNs and SNs on the fcth level of 
T* is denoted by . Let T be a feasible tree returned by the 
SRNP algorithm, and the set of RNs and SNs on the fcth level 
of T is denoted by . Thus, the ratio between the optimal 
solution and the solution returned by the SCA algorithm is 
given by 


u 






U +T-^) 

fc=l 

U +^.-5) 


( 15 ) 
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where I and I* denote the number of levels of T and T*, 
respectively. 

Because of the fact that a node can only be located on a 
certain level, i.e., P| = 0 and P| Lip, = 0 (i ^ j), 
equality ( fT5l l changes into 


U 

k^l 


U 


which leads to that 


E 

fc^l 

_ 

k^l 


(16) 




k^l 


< 


fc=l 






A:=l 


|L, 


+ 


I-^tI 


£ £ \^tA 

/c=l 


-/ 


E 


k 

rp. 

(17) 


Let OPTk be a minimum set cover for the fcth level of T, 
then, we can obtain that 


\Pp\ 


|L2| 


\Lir\ 


£ 


£ 14. 


£ 14- 


141 lo^p^il 

'^^^^'£|4.| 

i=l 


+ 


141 


fe=l 

IOPT 2 I 


l4i 


I OPT, I 


£ 14.1 

k=l 


£ 14.1 

fe=l 


(18) 


Since greedy algorithm is employed to solve the set covering 
problem that requires to select a subset of V{Lp~^) to cover 
the unconnected SNs Sk on the fcth level, considering that the 
approximation ratio of greedy set covering algorithm is In |4| 
El, we have that 


Vfc € {1, 2, 1}, 


141 

\OPTk\ 


<ln|4|. 


(19) 


In the first step of SCA, we check the existence of the 
feasible solution without RNs for DCRNP problem, thus, for 
the DCRNP problem input to the second step, its optimal 
solution contains at least one RN, i.e.. 


U ^T. ~ ^ 


fe=i 


> 1 , 


( 20 ) 


and according to equations (fTSI l- lfThl l this leads to that 


14 . ^ i5'i+1- 


fc=i 


( 21 ) 


Besides, due to the fact that \OPTk\ < 141 < |<S'|, 
combining with inequalities (fT9]l- (l2lT i. the formulation (fTsT i 


changes into 

141 \OPTi\ 

£ 14.1 

k=l 


141 


I0PT2I 


£ 14.1 

|5| + 1 
< Aln|S|. 


I4l \OPL\ 

£ 14.1 

fe=l 


( 22 ) 


Finally, the approximation ratio of SCA can be represented 
by 


U 4-5 


[j L^t^-S 


k^l 


0(ln|5|). 


(23) 


B. Time Complexity 

Let N = |5| + |C| + 1, then the time complexity of the 
shortest path tree algorithm is given by A^ Ig IIZTlI . The time 
complexity of the first step of the SRNP algorithm can be 
easily calculated, i.e., 0{N\gN). 

Then we analyze the time complexity of the Algorithmic In 
each iteration of the second step SCA, the main loop between 
lines 5-23 is executed. Therefore, we first analyze the time 
complexity to execute the main loop for once. In the inner loop 
between lines 9-14, the shortest path tree algorithm is applied 
to find the SNs fulfilling the inequalities (|9]) and dH. Thus, 
the time complexity of one iteration of this loop is 0{N Ig N). 
Additionally, this loop can be iterated at most N times, which 
leads to that the time complexity of this loop is 0{N'^ Ig A^)- 
Then, the greedy heuristic algorithm with a time complexity 
of 0(A(^) II 27 I is employed to find a minimum set cover for 
each level. Moreover, in each iteration of the second step of 
SCA, we will find the 1-hop neighbors within 0{N'^) running 
time. Hence, the time complexity of one iteration of the main 
loop is 

0{N^lgN)+0{N^) =0{N^lgN). (24) 

Since the main loop will be iterated at most A times, the time 
complexity of the second step of SCA is 0{N‘^ Ig A^)- 

In the main loop of the third step of SCA, the shortest path 
tree algorithm is implemented at most N times, and the time 
complexity to sort the remaining CDLs is O(A^). Thereby, 
the time complexity of the third step of SCA is 0{N^). Since 
the SCA algorithm is the combination of the three steps, the 
total time complexity, Tsca, of the SCA algorithm is given 
by 

TscA = 0{N\gN)+0{NHgN)+0{N^) 

= 0(Ar3), 

which indicates that the proposed SCA is a polynomial time 
algorithm. 


VII. Simulation Results 

In the simulations, SNs are randomly placed on a square 
field with side length of 100 units. To ensure a high probability 
that the DCRNP problem has a feasible solution, the amount 

































of CDLs is relative large, i.e., is set as 400, and these CDLs 
are randomly distributed on the deployment held. During the 
simulation, the amount of SNs is varying from 10 to 100. The 
simulation is carried out under both the homogeneous network, 
i.e., r = R = 10, r = R = 15 and r = R = 20, and the 
heterogeneous network, i.e., r = 10 and R = 15. SPTiRP 
algorithm proposed in m-ini is taken as the baseline. 
Since the heterogeneous network is not considered in mi- 
ED, when simulation is performed under the heterogeneous 
network, SPTiRP is adapted such that it can be implemented 
to the heterogeneous network. The notations used in the 
remaining part of this section are listed as follows: 

• iiiV(SPTiRP) and RN(SCA) represent the number of 
RNs deployed by SPTiRP and SCA, respectively. 

• i?A^(SPTiRP) — RN(SCA) represents the difference of 
deployed RNs between SPTiRP and SCA. 

A. Statistic Analysis 

The number of RNs deployed by SPTiRP and SCA under 
different conditions is shown in Fig. |6] and the difference of 
deployed RNs between SPTiRP and SCA (i.e., i?N(SPTiRP) — 
RN (SCA)) under different conditions is shown in Fig. |7] 

The results show that RNs deployed by SCA are fewer than 
those deployed by SPTiRP under different communication 
radii and delay constraints, which demonstrates the efficiency 
of SCA. It can be seen from Fig. |6(a)| that when the communi¬ 
cation radii of SNs and RNs is set as r = 10 and R — 15, SCA 
can save at most 5.4095 (5.4095/24.9524 ss 21.7%) deployed 
RNs comparing to SPTiRP. In Fig. |6(b)| the communication 
radii of SNs and RNs is set as r = i? = 20, SCA can 
save at most 0.8108 (0.8108/8.8108 Ri 10%) deployed RNs 
comparing to SPTiRP. In Fig. |6(c)| the communication radii of 
SNs and RNs is set as r = i? = 15, and the maximum saving 
due to SCA is 2.18 (2.18/17.5 sa 12.5%). In Fig. the 
communication radii of SNs and RNs is set as r = i? = 10, 
SCA can save at most 30.6667 (30.267/96.1429 Ra 31.48%) 
deployed RNs comparing to SPTiRP. 

It is clear in Fig. |6] that the number of RNs deployed 
by SPTiRP and SCA both increases with the decreasing of 
communication radii of SNs and RNs, and it is easy to 
explain this by that more RNs are required to connected the 
whole network as their communication radii become smaller. 
Fig. |6] also indicate that as the delay constraint is relaxed, 
SPTiRP and SCA will deployed fewer RNs under the same 
communication radii. This is due to that as the delay constraint 
is relaxed, the algorithm will select the solution containing 
fewer RNs, even though this will yield feasible paths with 
larger hop count for the given SNs. 

The number of RNs saved by SCA comparing to SPTiRP 
under different conditions is shown in Fig. Q It can be learned 
from these figures that as the communication radii becomes 
smaller, more RNs can be saved by SCA comparing to SPTiRP, 
i.e., the value of i?N(SPTiRP) - RN{SCA) grows. 

This phenomenon can be explained by that as the commu¬ 
nication radii becomes smaller, the topology becomes sparser 
and the hop count between the nodes on different paths 
becomes larger, therefore, the opportunity to merge different 


paths only by deleting the deployed RNs becomes smaller. 
Hence, the performance of SPTiRP deteriorates as the com¬ 
munication radii becomes smaller. On the other hand, SCA do 
not confront the limitation suffering by SPTiRP, therefore, the 
SCA algorithm will not suffer this deterioration. Consequently, 
the difference between these two algorithms grows larger as 
the communication radii becomes smaller. 

Additionally, we can see from Fig. [T] that as the delay 
constraint is tighten SCA can save more RNs than SPTiRP. 
This is due to that SCA can deploy fewer RNs on each level 
based on the set covering approach even when delay constraint 
becomes smaller. In contrast, when delay constraint is tighten, 
the RNs can be deleted from the original shortest path tree 
by SPTiRP become fewer, which reduce the performance of 
SPTiRP. As a result, the difference of deployed RNs between 
SPTiRP and SCA becomes larger. 

VIII. Conclusion 

In this paper, the DCRNP problem is studied. For the 
NP-hard nature of the DCRNP problem, an approximation 
algorithm-SCA is proposed to handle this problem based on 
the set covering method. The novelty of the proposed SCA 
algorithm lies in its multi iteration procedure, and in each 
iteration the set covering method is applied to select the SNs 
or CDLs. The approximation ratio of the SCA algorithm is 
proved to be ©(Inis'!). The time complexity of the SCA 
algorithm is also rigorously analyzed, and it shows that the 
SCA is a polynomial time algorithm with the time complexity 
of 0{N^). Finally, extensive simulations are carried out to 
evaluate the performance of the SCA algorithm. The simula¬ 
tion results show that the SCA algorithm can significantly save 
the deployed relay nodes comparing to the existing algorithm 
within a worthy running time. 
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